Forecasted 21 cm constraints on compensated isocurvature perturbations 
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A "compensated" isocurvature perturbation consists of an overdensity (or underdensity) in the 
cold dark matter which is completely cancelled out by a corresponding underdensity (or overdensity) 
in the baryons. Such a configuration may be generated by a curvaton model of inflation if the cold 
dark matter is created before curvaton decay and the baryon number is created by the curvaton decay 
(or vice- versa). Compensated isocurvature perturbations, at the level producible by the curvaton 
model, have no observable effect on cosmic microwave background anisotropies or on galaxy surveys. 
They can be detected through their effect on the distribution of neutral hydrogen between redshifts 
30 to 300 using 21 cm absorption observations. However, to obtain a good signal to noise ratio, very 
large observing arrays are needed. We estimate that a fast Fourier transform telescope would need 
a total collecting area of about 20 square kilometers to detect a curvaton generated compensated 
isocurvature perturbation at more than 5 sigma significance. 



I. INTRODUCTION 

Current observations can be fit by a simple 6 param- 
eter model of a flat ACDM Universe with a nearly scale 
invariant, adiabatic, and Gaussian primordial perturba- 
tion spectrum (sec for example [HQ). The primordial 
perturbations may have been generated by a period of 
accelerated expansion in the early Universe, known as 
inflation (see for example (H). Inflation can be driven by 
a potential dominated scalar field (the inflaton) . The cur- 
rent Universe contains cold dark matter (CDM), baryons 
(£>), neutrinos (y) and photons (7). They can arise from 
the decay products of the inflaton. We assume the recent 
accelerated expansion of the Universe is been driven by 
a cosmological constant which is unperturbed. Density 
perturbations in the other components can be inherited 
from the density perturbations in the inflaton which arise 
from vacuum fluctuations that are amplified to scales 
larger than the Hubble horizon during inflation. The 
perturbations are known as adiabatic if 
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where Sp is the perturbation in the density, a dot indi- 
cates a derivative with respect to time and i and j are 
each one of (CDM, b, v, and 7). Non-adiabatic (also 
known as isocurvature or entropy) perturbations cannot 
arise if all the constituents of the Universe are the result 
of the decay of a single inflaton (see for example 0, [f| ) . In 
order to generate isocurvature perturbations there has to 
be more than one light degree of freedom present during 
inflation. In the curvaton model [6, 7] , the perturbations 
generated by the inflaton are negligibly small. But, there 
is a second light field (known as the curvaton) which also 



acquires perturbations, from vacuum fluctuations, dur- 
ing inflation. At the end of inflation, the inflaton decays 
into radiation whose density decreases like a -4 , where 
a = 1/(1 + z) is the scale factor and z is the redshift. 
After inflation, when the Hubble parameter drops below 
the curvaton mass, the curvaton oscillates in the well of 
its potential and its density decreases like a~ 3 and so its 
energy density relative to the radiation increases as a. 
When the Hubble parameter drops below the curvaton 
decay rate, the curvaton decays and its decay products 
are responsible for the observed density perturbations. 
If all of the constituents of the current Universe origi- 
nate from the curvaton, then the resulting perturbations 
will be adiabatic. However, if some of the constituents 
originate from before the curvaton has non-negligible en- 
ergy density and some from the curvaton, or its decay 
products, then isocurvature modes are possible Q. The 
magnitude of the adiabatic mode can be given in terms of 
the curvature perturbation in the constant density gauge, 
£ (see for example [!, [ij [l(| ) ■ In the cases we consider 
it is virtually equivalent to the curvature perturbation 
(or its negative, depending on the sign convention) in 
the comoving gauge, TZ. The isocurvature term between 
constituent i and j can be quantified by 
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where H is the Hubble parameter, the second equality 
follows from the continuity equation, and w is zero for 
baryons and CDM and a third for photons and neutri- 
nos. If the CDM is created just after inflation, before 
the curvaton has non-n egligibl e energy density (as in the 
case of Wimpzillas fill [l2l.Tl3j). then Q 
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Such a large negatively correlated isocurvature perturba- 
tion on its own is ruled out by current data However, 
it may be offset by a positively correlated isocurvature 



2 



perturbation between the baryons and the photons. If 
the curvaton decay generates the baryon number (e.g. 
baryon number from the decay of a right-handed sneu- 
trino curvaton [15j). then Q 
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where r f=a Pcurvaton/ptotai at the time of curvaton decay. 
Observations of the anisotropy in the cosmic microwave 
background (CMB) are unable to distinguish between a 
Sb,-y and a Scdm^ isocurvature perturbation [l4j]. Also, 
as discussed in Sec. [HI they cannot be distinguished with 
galaxy redshift surveys either. However, the effect of a 
CDM isocurvature mode is f2cDivi/^b times larger than 
that of a baryon isocurvature mode with the same ampli- 
tude. Where Qi is the density of component i divided by 
the critical density. It follows that there is a, currently 
unconstrained, compensated mode where 
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The case where the CDM is created before curvaton de- 
cay and the baryon number by curvaton decay corre- 
sponds to 1 = —3. 

As pointed out in [JJI, [l7|, high redshift 21 cm obser- 
vations can distinguish between Scdm,7 and Sb. y In this 
article we estimate how large an array would be needed 
to detect a curvaton generated compensated isocurvature 
mode. The layout of the article is as follows. In Sec. [TTT 
IIII1 and IIV1 the effect of compensated isocurvature per- 
turbations is discussed. In Sec. |V] forecasts for how well 
the compensated isocurvature perturbations can be con- 
strained are made for possible future 21 cm experiments. 
The conclusions and discussion are given in Sec. IVII 



II. EFFECT OF COMPENSATED 
ISOCURVATURE ON THE CMB 

In [3] it was shown that the compensated mode is not 
observable in the CMB anisotropy. Here we go through 
the argument in a bit more detail and highlight the ef- 
fect of baryon pressure. The temperature perturbation 
in the CMB is given by the sum of the temperature per- 
turbation from the adiabatic mode and that from the 
compensated mode: 
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The compensated isocurvature mode has all perturba- 
tions set to zero except the baryon and CDM densities 
which have SpcuM — —&Pb and so the total density is un- 
perturbed. In the compensated mode, the divergence of 
the baryon fluid velocity is sourced by the baryon density 
(see for example [U[). But, the sourcing term is of the 
form c 2 s k 2 5pb/ Pb where the baryon sound speed is give by 
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and Tb is the baryon temperature (which is equal to the 
photon temperature until z ~ 300) and p is the mean 
molecular weight (including free electrons and all ions 
in H and He). Only wave numbers for which k 2 c 2 s > Ti 2 
(where Ti is the comoving Hubble parameter) are affected 
by the baryon sound speed. At last scattering, this corre- 
sponds to k > 200/iMpc -1 which corresponds to angular 
scales of the observed CMB of £ > 2 x 10 6 which is many 
orders of magnitude greater than the observable CMB 
intrinsic anisotropy which due to a combination of Silk 
damping and foregrounds is unobservable for £ > 3000. 
Note that c 2 is not the square of the photon/baryon fluid 
sound speed which is of order a 1/3 during tight cou- 
pling, where we are using units with the speed of light 
set to unity. Before last scattering, the divergence of 
the photon fluid velocity is sourced by the divergence of 
the baryon velocity. The photon's density and velocity 
perturbations are not sourced by the baryon density per- 
turbations. Hence, except on negligibly small scales, the 
compensated isocurvature mode has all perturbations re- 
maining zero except for Spb and SpcuM which are initially 
non-zero but do not grow as over all there is no pertur- 
bation to the metric. It follows that the observed CMB 
fluctuations are unaffected by the presence of a curvaton 
generated compensated isocurvature mode. 



III. EFFECT OF COMPENSATED 
ISOCURVATURE ON GALAXY SURVEYS 

The dimensionless power spectra is given by 

V = (k 3 /2n 2 ) (Sp(k) 2 ) Ip 2 . (9) 

It is plotted for the baryon and CDM perturbations in 
Fig. [1] for scales smaller than the Hubble horizon. The 
power spectra for this and the other figures in this ar- 
ticle were evaluated using the CAMB sources [43[ pro- 
gram which solves the Boltzmann equation for the pho- 
ton distribution function numerically [T^ |. Our fiducial 
adiabatic model has WMAP5 maximum likelihood pa- 
rameters [TJ]: tt b h 2 = 0.0227, n CDM h 2 = 0.108, n = 
0.961, r = 0.089, P C \ k= ^ hMp c^ = 2Al * ^~ 9 ,h = 
0.724 where H\ z=0 = 100/ikm-s _1 -Mpc _1 , n is the spec- 
tral index, and r is the optical depth. The compensated 
mode for the curvaton model is taken to be T = —3 
which corresponds to the case of the CDM being created 
before curvaton decay and baryon number being created 
by curvaton decay. The compensated mode where baryon 
number is created before curvaton decay and CDM by 
curvaton decay corresponds to X = 3f2b/O c = 0.63 and 
is harder to detect. 

The compensated mode density perturbations are de- 
termined by the isocurvature perturbation, Eqs. ^ and 



3 




k/(hMpc _1 ) 



10" 


1 


10" 


-2 


10" 


-3 


10" 


-4 


10" 


-5 


10" 


-6 



10~ 3 10~ 2 10 _1 10° 10 1 10 2 
k/ChMpc" 1 ) 




k/(hMpc _1 ) 

FIG. 1: The evolution of the dimensionless power spectrum 
of the baryons (red) and the CDM (blue). The adiabatic 
(solid) and adiabatic + curvaton compensated isocurvature 
mode (dashed) are also plotted. The baryon and CDM power 
spectrum have been scaled by fi(,/fi m and ficDM/^m respec- 
tively, where fi m = Qcdm + fit- 



(|rj|), which in the curvaton model has the same primor- 
dial spectral index as the curvature perturbation (V oc 
k n ~ 1 ). The compensated mode density perturbations do 
not change with time except on scales smaller than the 
baryon pressure scale corresponding to k ~ 200/iMpc -1 
at the rcdshifts probed in Fig. [T] For scales larger 
than the Hubble horizon at radiation/matter equality 
(k < O.Ol/iMpc -1 ) the adiabatic mode density perturba- 
tions are determined by the curvature perturbation and 
from Poisson's equation have the form V oc k n+3 / (1+z) 2 . 
For larger k, there is a reduction in the slope due to 
the radiation pressure suppressed growth of sub-horizon 
scales during the radiation era. During the matter era 
all scales grow as (1 + z) -1 in the adiabatic mode den- 
sity perturbation. It follows that the contribution of the 
compensated mode to the square root of the normal- 
ized power spectrum decreases with scale and redshift 



as k~ 2 (l + z) for k < 0.01/iMpc -1 and similarly but with 
slightly less negative slope in k for k > 0.01/iMpc -1 . 
The total density perturbation is the sum of the com- 
pensated mode and the adiabatic mode density pertur- 
bations. The CDM compensated mode density pertur- 
bation is completely negatively correlated with the adi- 
abatic mode CDM density perturbation on large scales. 
The downward spike seen in the y/V between k = 10~ 3 
and 10~ 2 ft.Mpc is where the adiabatic and compen- 
sated mode CDM density perturbations completely can- 
cel each other. 

The number of gravitationally collapsed halos of 
mass M is determined by the variance of the non- 
relativistic matter field (evaluated using linear the- 
ory) smoothed with a top hat filter on a scale R ~ 

/i- 1 Mpc(M/lO 12 M ) 1/3 [H. The number of objects 
is exponentially suppressed when the corresponding 
smoothed variance is less than about 1.686 2 [20]. The 
variance of the smoothed matter field can be approxi- 
mately read off from the dimensionless power spectrum 
at scale k ~ 1/R. The first stars form in halos of mass 
10 6 M Q . This corresponds to k ~ 10 2 /i _1 Mpc On the 
other end of the scale, clusters of galaxies can have masses 
up to about 10 15 Af Q which correspond k ~ 10/i _1 Mpc. 
and tend to form about z ~ to 2. As can be seen from 
Fig. [TJ the effect of the compensated mode is negligible 
for k > 0.1ft _1 Mpc for z < 100. It follows that a curva- 
ton generated compensated isocurvature mode will have 
a negligible effect on collapsed gravitational structures 
and so is not detectable by galaxy surveys either. 

IV. EFFECT OF COMPENSATED 
ISOCURVATURE ON 21 CM OBSERVATIONS 

Absorptions in the observed CMB at wavelengths 

A = 21.1(1 + z)cm (10) 

can be used to probe the distribution of neutral hydrogen 
in the early Universe (see for example plT ] for a review). 
The amount of absorption of the CMB is maximized at 
z ~ 70 where there is the biggest difference between the 
CMB temperature and the spin temperature of the neu- 
tral hydrogen. If there is an overdensity of baryons, and 
hence neutral hydrogen, then there will be more absorp- 
tion. The amount of neutral hydrogen also affects how 
the spin temperature evolves with time. There are many 
other subtle effects, such as distortions due to the pe- 
culiar velocities of the neutral hydrogen, but these will 
be subdominant at the bandwidths, rcdshifts and length 
scales we look at here [l?], [HJ . To a first approximation, 
the 21 cm observations can be thought of as mapping 
the density of neutral hydrogen (or effectively baryons) 
for redshifts in the range 30 < z < 300. For lower red- 
shifts, 21 cm can be used for probing reionization (see for 
example (23j ) and as a means of doing galaxy surveys (see 
for example [lOj]). Also, for z < 30, the neutral hydro- 
gen can still be mapped to constrain the matter power 
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FIG. 2: The 21 cm angular power spectrum for the adiabatic 
mode is plotted (black, solid) with one sigma error bars, using 
At — 1. The signal from the adiabatic plus the curvaton 
compensated isocurvature mode is also plotted (dashed, red). 
The error bars are for a 20 km 2 FFT telescope and a band 
width of 8 MHz. They start at the minimum value of £ which 
is detectable once foregrounds have been taken into account. 



spectrum, provided one marginalizes over the reioniza- 
tion model (see for example [251 HH, These are the 

goals of current (e.g. MWA J4I, LOFAR [H]) and the 
next generation (e.g. SKA j46[) surveys. But next+1 
generation surveys may probe the dark ages for z > 30. 
As the neutral hydrogen follows the baryon distribution 
at high redshifts it will have a signal for a compensated 
isocurvature mode 0, In Fig. [5] the 21 cm signal 
is plotted at different redshifts for the adiabatic and the 
adiabatic plus compensated isocurvature mode. The er- 
ror bars are for those scales measurable after foregrounds 
have been removed and will be discussed in Sec. [Vj For 
£5r/ r ^> 1 , (corresponding to £ ^> 20 in Fig. \2§ the an- 



gular power spectrum can be approximated by [26 

l(l+l)Ct/27T<XjV b \ k=i/r 



(11) 



where r is the comoving distance to the center of the sur- 
vey, Sr is the comoving width of the survey, and Vb is the 
dimensionlcss power spectrum of the baryons. The extra 
factor of l/£ accounts for the smoothing effect of the sur- 
vey window. For the redshifts in Fig. [21 r ss 10 4 h Mpc -1 . 
Fig. [2] shows that the signal to noise is greatest where 
£ ~ 300, which corresponds to k ~ 0.03 ft. Mpc -1 . As can 
be seen in Fig. [1] the baryons in the compensated mode 
add non-negligible power at these scales and this can be 
seen in the corresponding compensated isocurvature plus 
adiabatic mode curves in Fig. [5] 



V. FORECASTS 

The FFT telescope [29[ is well suited for measuring 
high redshift 21 cm. It has antennas arranged in a regu- 
lar grid allowing an FFT to be used when calculating the 
correlations between antennae. This greatly reduces the 
associated computational cost, which otherwise becomes 
prohibitive for large arrays with many antennae. Addi- 
tionally, its large field of view means it can be used in 
drift mode allowing a quarter of the sky to be surveyed. 
As seen in Fig. 8 of [20|, for the large collecting area 
needed to survey a redshift of z = 50, the FFT telescope 
is the cheapest option. The noise power spectrum is [29( 
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C noise d- 
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(12) 



(13) 



(14) 



(15) 



where the beam function is taken to be Gaussian 

Bf = e° 2e2 

with the resolution given by [29| 
6 = X/VA 

where A is the FFT telescope area. Also, [29j 

^rnoise _ 47T A_Aky^ys 
" v jcover ^Q cr ' 

Here / s k y = il ma p/47r is the fraction of the sky covered by 
the map, f2 is the field of view, we have introduced the di- 
mensionless parameter rj = Au/u = Av c/A to denote the 
relative frequency bandwidth, and c is the speed of light. 
In CAMB sources the averaging over frequency is done 
with a Gaussian of standard deviation Ai//(2y / 7r) [3(| 
specified by "rcdshift_sigma_Mhz" in the CAMB sources 
initialization file. The observation time is denoted by r, 
the system temperature by T sys , and J covor is the fraction 
of the area covered by the array antennas. 

The one sigma error bars for the angular power spec- 
trum are [291 



AC 1 



{21 + l)A£f sky 



(Ci + C\ 



(16) 



20 



FIG. 3: The Fisher information of the compensated isocurva- 
ture parameter as a function of £ is plotted for a FFT telescope 
with an area of 20 km 2 . Redshift 30 (solid, blue), 40 (dashed 
green), and 50 (dotted red) are shown. 



Removing foregrounds at these high redshifts will be ex- 
traordinarily challenging but, if the foregrounds are suf- 
ficiently smooth in frequency, it may well be possible. 
We will assume that foregrounds can only be removed 
for t > train where ^ m j n = fc m j n r with the minimum wave 
number corresponding to the bandwidth and given by 

m 

fc min = 2ir/Ar . (17) 
The Fisher information matrix is given by [32j 



where 



F, 



1 dC e dd 
AC 2 d Pl dp 3 



(18) 



(19) 



and pi is parameter i. 

We set the observation time to r = 365 x 24 hours and 
following [33[, T sys = 200K x [(1 + z)/10] 2 - 6 , the field of 
view to be O = O map = tt, and / covcr = 1. We set the 
bandwidth to be Av = 8 MHz. The Fisher information 
for the compensated isocurvature parameter is plotted 
as a function of I in Fig. [31 As can be seen, at low 
£, redshift 40 measurements are more sensitive while at 
higher I redshift 30 measurements are more sensitive. 

In order to remove degeneracies with the other parame- 
ters, we include PLANCK [13] CMB forecasts. These do 
not constrain the compensated isocurvature parameter 
(I), but they do help to constrain the other parameters 
which may be degenerate with I. We use the forecasted 
temperature and polarization measurements from the 70, 
100,143 and 217 GHz bands of PLANCK, see Sec. Ill 
of [34[ for more details. The combined Fisher matrix 
(-p2icm+PLANCK)) is obtained by adding the 21cm and 
PLANCK Fisher matrices. The forecasted covariancc 
matrix for the parameters is given by inverting the com- 
bined Fisher matrix. The number of "sigma" for which 
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FIG. 4: The number of sigma that the curvaton compensated 
CDM isocurvature mode would be detected is plotted. Each 
point in the plot is for a 8 Mhz band width experiment. 



a curvaton compensated isocurvature mode could be de- 
tected is given by 3/ (^"lcm+PLANCK) )zz where subscript 
II denotes the element in the row and column corre- 
sponding to the compensated isocurvature parameter. A 
contour plot of the number of sigma detection, as a func- 
tion of the area of the FFT telescope and the redshift 
probed, is plotted in Fig. 0J As can be seen, in order 
to obtain a five sigma detection, a redshift around 40 is 
optimal and a FFT telescope with an area of at about 
20 km 2 would be needed. According to [29], that area 
FFT telescope would cost of order several billion US dol- 
lars. However the need to place such an array beyond the 
Earth's ionosphere, for example on the moon [35j, would 
raise this cost significantly. 

The parameter with the most degeneracy with the 
compensated isocurvature mode is the baryon density. 
As can be seen from Fig. [5j there is a difference in 
the slope of the degeneracy for redshifts 30 compared 
to higher redshifts. This may indicate that combining 
redshifts will improve the constraints but the results will 
probably not be significantly improved as the effective 
area of the telescope can only be optimized for a partic- 
ular redshift and will then degrade for other redshifts. 



VI. CONCLUSIONS 

We have shown that it is in principle possible to use 
21 cm measurements to detect compensated isocurvature 
perturbations produced by the curvaton. At redshifts of 
about 40 the baryon perturbations are sufficiently differ- 
ent from the adiabatic case to give a detectable 21 cm 
signal, at the 5 sigma level, provided a FFT telescope of 
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FIG. 5: The contours containing 95% of the probability distri- 
bution for a 20 km 2 FFT telescope combined with a PLANCK 
CMB experiment. The axes are for the CDM compensated 
isocurvature mode and the baryon density. All other param- 
eters are marginalized. Redshifts 30 (solid blue), 40 (dashed 
green), and 50 (dotted red) are plotted. The baryon density 
has the greatest degeneracy with the compensated isocurva- 
ture mode. 

about 20 km 2 or larger is used. 

The curvaton model can also produce non-Gaussianity 
of magnitude f n i = 5/4r. So for the compensated isocur- 
vature value of r f=a 0.17 a f u \ « 7.35 will be generated. 
This may be detectable by the PLANCK satellite [HI 
and potentially even 21 cm observations [37| . 

In this article we have looked at a curvaton generated 
compensated isocurvature mode which from using Sp 7 = 
in Eq. ([2]) and substituting into Eq. ((6|) gives 

^ ~ 3C(0 CDM /0 6 ) ~ 10- 3 
Pb 

on large scales. If compensated modes which have 
3pb/pb ~ lO^ 1 arc considered, then it becomes pos- 
sible to use many more observational constraints such 
as the scatter in light element abundances due to in- 
homogeneous big bang nuclear synthesis, galaxy cluster 
gas fraction measurements, B-mode polarization induced 



from inhomogencous rcionization, and possibly even very 
large galaxy surveys or high redshift QSO counts [H| . It 
would be interesting to try and find early Universe mod- 
els that could induce such a large compensated isocurva- 
ture mode. 

In the case of the curvaton generated compensated 
isocurvature mode, there is also likely to be some resid- 
ual correlated CDM (or equivalently) baryon isocurva- 
ture modes unless there is some specific reason for the 
baryon and CDM isocurvature modes to exactly cancel 
each other out as in the condition specified in Eqs. (O 
and ©. The WMAP3 data 95% confidence constraint 
was [39| —0.42 > Sb,j/C < 0.25. Curvaton model con- 
straints have also been evaluated in WMAP5 @,H3], but 
they have effectively assumed Sb,-y > and so we can't 
directly use them. If the compensation in Eq. ([5]) is not 
exact, there will be a residual 

^ 1 = _3^cdm +3 1_t_ 
C *h r 

Using the WMAP3 constraints in the above formula we 
find that r would need to be within about 4% of the com- 
pensated value of r = Ofe/OcDM ~ 0.17. The updated 
constraints from WMAP5 are likely to be even stronger. 
So the fine tuning will probably be more severe than 4%. 
However, it is intriguing that the baryon density is so 
similar to the CDM density. This may be a hint that 
baryogenesis and CDM creation are in some way related 
and some models have been proposed to account for this 
(see for example [U H^j). It would be interesting to 
investigate if there were some variant of these models 
that could naturally produce a sufficiently compensated 
isocurvature perturbation to be consistent with current 
constraints. 
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